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1 Introduction 


Constants of motion (invariants or first integrals) are of central importance 
to study dynamical systems in general. This has motivated the development 
of several procedures by which constants of motion for dynamical systems 
(see [1] for a recent review) can be constructed. Among these, the appli¬ 
cation of Noether’s theorem [2, 3] deserves special attention in view of the 
physical appealing of the method. Indeed, Noether’s theorem provides links 
between continuous symmetries of the action functional and constants of mo¬ 
tion of the system. Among the classical results obtained by application of 
Noether’s theorem are the conservation of energy, linear momentum and an¬ 
gular momentum. The associated symmetries correspond to the invariance 
of the action functional under time translation, space translation and space 
rotations, respectively. 

Recently, Ermakov systems [4] were analyzed [5]—[8] via Noether’s theo¬ 
rem. Ermakov systems, which are pairs of coupled, second-order ordinary 
differential equations, always admit at least one constant of motion, the 
Ermakov invariant. The work developed in [5]-[8] provided a consistent ex¬ 
planation for both the equations of motion and the invariant of Ermakov 
systems, as a result of an underlying Noether symmetry. The approach 
proposed, however, applied only to Ermakov systems where one of the equa¬ 
tions decoupled from the other. This limitation has motivated comments 
[9] concerning the difficulty of deriving the Ermakov invariant from symme¬ 
try considerations only. A closer examination shows that references [5]-[8] 
assumed that one of the variables in the Ermakov system plays the central 
role, the other being regarded as a mere auxiliary variable. This point of 
view naturally led to systems where one of the equations (the so called aux¬ 
iliary equation) decoupled from the other, of a more important character. In 
the present work, we consider the Ermakov systems as essentially coupled, 
with no dynamical variable playing any distinguished role. For this, we pro¬ 
pose for the Ermakov systems not a one-dimensional but a two-dimensional 
Lagrangian formalism. Using this Lagrangian description, we search to con¬ 
struct coupled Ermakov systems for which Noether’s theorem applies. In this 
context, both the Ermakov and all additional invariants - when they exist - 
should result from the associated Noether symmetry which then provides a 
general framework to define, classify and solve Ermakov systems. We also 
remark that the use of a two-dimensional Lagrangian formalism for Ermakov 
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systems has been suggested in a recent work [10]. 

A central issue with Ermakov systems concerns their complete integra- 
bility. Specihcaly, we analyse the complete integrability of Ermakov systems 
in the Liouville sense [11], which, in turn, requires the existence of an under¬ 
lying Hamiltonian structure. In fact, complete integrability requires, in this 
cases, the existence of at least a second invariant. The existence of Ermakov 
systems with a second invariant was first addressed by Goedert [12], who 
obtained a class of Ermakov systems possessing a second constant of motion. 
For Hamiltonian Ermakov systems, Liouville’s theorem [11] warrants com¬ 
plete integrability if a second constant of motion exists in involution with the 
Ermakov invariant. Notice also that, strictly speaking, Liouville’s theorem 
requires compact level surfaces of the constants of motion. 

The present work focus on the application of Noether’s theorem to Er¬ 
makov systems. The Lagrangian formalism constitutes the usual framework 
for the application of Noether’s theorem [2] (and references therein). We 
thus need Ermakov systems that may be derived from a Lagrangian formal¬ 
ism. This has already been dealt with in the literature [13, 14]. Cervero 
and Lejarreta [13] showed that conventional Hamiltonian Ermakov systems 
are completely integrable. By conventional Ermakov systems, we mean Er¬ 
makov systems having a frequency function depending on time only. How¬ 
ever, as pointed out previously [15], a frequency function depending also on 
the dynamical variables does not destroy the main feature of Ermakov sys¬ 
tems, namely the existence of the Ermakov invariant. This has motivated a 
number of studies [12], [14, 16, 17] on Ermakov systems with a generalized 
frequency function, depending on dynamical variables besides time. Follow¬ 
ing this trend, we do not restrict the frequency functions dependence on 
time only. Nevertheless we will show that the integrability of conventional 
Hamiltonian Ermakov systems results from a Noether symmetry, a fact that 
was not clear at all in the original reference [13]. Finally, we remark that 
dynamical symmetries have been considered in reference [18] for Ermakov 
systems of conventional type and not necessarily having Lagrangian charac¬ 
ter. Let us stress that in the present work the interest is centered on Noether 
symmetries, which may be assessed only when a variational description of 
the system is available. Besides, we do not restrict ourselves to conventional 
Ermakov systems. 

The paper is organized as follows. In section 2, we obtain a class of 
Ermakov systems admitting a Lagrangian description, directly from the class 
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of Hamiltonian Ermakov systems derived in reference [14]. In section 3 we 
show, by use of the converse of Noether’s theorem [2], that the Ermakov 
invariant is the result of a Noether dynamical symmetry. Within the class 
of Lagrangian Ermakov systems, we also find the subclass of systems that 
possess Noether point symmetry. The Lagrangian Ermakov systems with 
Noether point symmetry, having a second constant of motion in involution 
with the Ermakov invariant, are Liouville integrable. 

2 Lagrangian Ermakov systems 

Ermakov systems are generally [14] defined as the system of equations 


X + u 2 x = 

-^2 > 
yx z 

( 1 ) 

y + uj 2 y = 

-\g(x/y) , 

xy z 

( 2 ) 


where / and g are arbitrary functions of the indicated arguments and uj is an 
arbitrary function of time, the dynamical variables x and y, and their time 
derivatives. For physical reasons, it is convenient to restrict the dependence 
to u> = u(t,x,y,x,y). For brevity, we call u the frequency function, even 
if it is not associated with any characteristic frequency of the system. We 
recall that, as shown in [14], only one arbitrary function of y/x is sufficient 
to specify a generic Ermakov system if the frequency function is allowed to 
depend on the dynamical variables. However, for easier comparison with the 
results in the literature, we take (1-2) as our definition of Ermakov systems. 

The central property of Ermakov systems is their Ermakov invariant, 

1 = \( X V - yx ) 2 + J f(X)dX + J g(X)dX. (3) 

The Ermakov invariant is constant along the trajectories of (1-2) indepen¬ 
dently of the detailed dependence of c o on its arguments. The Ermakov 
invariant is most useful to construct a nonlinear superposition law [15] relat¬ 
ing the solutions of (1) and (2). Here, however, we are interested in Ermakov 
invariants as one ingredient for the derivation of completely integrable Er¬ 
makov systems. 
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Given the variational formulation, we can use some more specialized tools 
to investigate integrability, such as Noether’s theorem [2]. This fact makes it 
interesting to consider not the general Ermakov systems (1-2), but only their 
subclasses that admit a variational formulation. In particular, we shall focus 
on Lagrangian descriptions, since it is on configuration space that Noether 
theorem is usually formulated [2], 

We, therefore, pose the following question: under what conditions on the 
functions u>, f and g does (1-2) possess a variational description? This ques¬ 
tion was partially answered by Cervero and Lejarreta [13], who considered 
Ermakov systems with usual frequency functions, and by Haas and Goedert 
[14], in the case of frequency functions depending on dynamical variables. 
Both works adopted a variational description of the Hamiltonian type. Here, 
however, we are specially interested in Lagrangian descriptions and their 
connection with Noether’s theorem. Therefore, we rewrite the results of ref¬ 
erence [14] in terms of a Lagrangian formalism. Here this is the more direct 
approach, even if the usual procedure is to derive the Hamiltonian formalism 
from the Lagrangian description. 

Using the results of [14] we may consider (1-2) as obtained from a La¬ 
grangian function that represents the sum of a quadratic kinetic energy term 
and a potential energy term provided that the frequency function is taken in 
the form 


u> 


2 


\dV_ 

RdR 


(.R, t ) + 


R A ’ 


(4) 


where a (9) is defined as 


a(6) 


(.A cos 9 + B sin 6) 
i/j 2 (9 ) sin 9 cos 2 9 


/(tan#) + 


(B cos 9 + C sin 9) 
i/j 2 (9 ) sin 2 9 cos 9 


9 (cot 9) 



(5) 


/ and g being kept totally arbitrary. In (4) V(R, t ) is an arbitrary function, 
and, R and 9 are defined by 

R 2 = Ax 2 + 2Bxy + Cy 2 , 9 — arctan(y/x), (6) 


with A, B, C and k numerical constants satisfying 

k = AC - B 2 ± 0 . 


(7) 
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Also, in (4) we used the notation 


ip 2 [O') = (a cos 2 9 T ‘IB sin 6 cos 9 -\- C sin 2 9^j . (8) 

As mentioned before, in general to does not represent any physical frequency 
of the system, but we keep the nomenclature for convenience of notation. 
The associated Lagrangian function is then 

L = T-V , (9) 


where the kinetic energy is 

T = ^(Ax 2 + 2Bxy + Cy 2 ), (10) 

and 

— K ( /“tan 6 rcot 9 \ 

V = V(R, t) + ^U f ^ dX + J sWdXj . (11) 

For the sake of completeness, we writte the corresponding Hamiltonian, 

H = — (Cp^ + 2 Bp x p y + Ap 2 ) + V . (12) 

An important point concerning Ermakov systems is their linearisation. 
As shown in [19], Ermakov systems with frequencies depending only on time 
are always reducible to a linear second-order ordinary differential equation. 
Latter on [20], the same linearisation was shown to apply to a subclass of Er¬ 
makov systems with frequency functions depending on dynamical variables. 
In the final part of the section, we verify in what cases the Lagrangian Er¬ 
makov systems just described may be linearisable following the procedure 
of references [19, 20]. This procedure relies on the introduction of a new 
dependent variable p given by 


c p = a(t)/R , (13) 

where a(t) is a function of time to be choosen conveniently. Also, we must use 
the polar angle 9 as the new independent variable. The change of independent 
variable is acomplished by the relation 

6 = h(9-I)/B 2 , (14) 
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( 15 ) 


which follows from the Ermakov invariant, where 

\/o ( /'tan# /'cot# \ 

h (<?; i) = -^{i-j + J g ^ dX ) ■ 

Actually, in references [19, 20] the variables used are a(t)/r and 0, where 
(r, 0) are the conventional polar coordinates, but here the choice (13-15) can 
be shown to be more effective. 

Using the Lagrangian equations of motion obtained from (9) and the 
relations (13-15), we get 


2 dV , , dhdl P to t _ 2 

h Hm +h mHi) +2Klv ~ a 




(16) 


where y is U expressed in terms of the new dependent variable, 


X = x{v,t) = V(a/ifi,t). 


(17) 


The left-hand side of eq. (16) is a linear form. Obviously, the complete 
equation will have a linear character if and only if the right-hand side of (16) 
is a linear form, which imposes 


X = xif) + 



bip 


act 

V ’ 


(18) 


where y(f) is an arbitrary function of time which can be set to zero with¬ 
out loss of generality since it has no influence on the equations of motion. 
Moreover, a, b are arbitrary numerical constants. The resulting linearisation 
reads 


i2 jV 
dd' 2 


+h ^ +2[Kl+a)v 


b. 


(19) 


The reconstruction of the solution of the original Ermakov system from the 
solution of (19) is made following the steps described in [19]. 

Here, the important point is that not all Lagrangian Ermakov systems 
arc linearisable in terms of the proposed change of variables. Indeed, (18) im¬ 
poses a particular dependence of the originally arbitrary function V, namely 

- . „ . a b a R 2 . . 

(20) 

For a = cte., the resulting equations of motion can be shown to be in the 
class of Kepler-Ermakov systems [21]. For general a(t), we have a generalized 
time-dependent Kepler problem [22], 
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3 Noether symmetries 

Let us consider the transformations 


x = x + erji (x,y,x,y,t), 

y = y + ey 2 (x,y,x,y,t ), (21) 

t = t + er{x,y,x,y,t), 

where £ is an infinitesimal parameter and r/i, r] 2 and r are functions to be 
determined. As these functions may depend on the velocities, we are dealing 
with transformations of dynamical type. Dynamical symmetries may also 
contain a dependence of r, r/i or r/ 2 on higher derivatives, but, for simplicity, 
we restrict the treatment to symmetries of form (21). 

To proceed we introduce the following more convenient notation, 


q ={x,y), V = {VuV2):i 

and consider the action functional 

S[q] = [ Ldt. 

Jto 


( 22 ) 


(23) 


Noether’s theorem establishes a correspondence between continuous trans¬ 
formations that leave the action functional invariant (up to an additive nu¬ 
merical constant) and conservation laws. In particular, Noether’s symmetry 
criterion [2] states that 


d L d L dL 

r aF + "'^ + ( "^ Tq) 'a5 


r L = A , 


(24) 


for A = A(q, q, t) a function of the indicated arguments. The associated 
Noether invariant is 


J = T 



(25) 


The Noether condition (24) may be used to determine simultaneously 
both the invariant and the action symmetry. However, if one invariant is 
known, we can always obtain an associated Noether symmetry, by use of the 
converse of Noether’s theorem [2], This usually leads to a symmetry of a 



dynamical character. In the present case, the Ermakov invariant is available, 
and the procedure may help to understand precisely its symmetry origin. 
Hence, we briefly review the converse of Noether’s theorem, and apply it to 
our Lagrangian Ermakov system. 

The converse of Noether’s theorem [2], states that if /(q, q, t) is an invari¬ 
ant for a Lagrangian dynamical system with regular Lagrangian L(q, q, t), 
then it has the Noether symmetry (using component notation and summation 
convention) 

i = q‘ + r-(-g't^ + Tq^ , (26) 

t = t + £T, (27) 

where r = r(q, q, t) is an arbitrary function and g lJ is the inverse of the 
Hessian matrix of the Lagrangian, 

msj 1 - <28> 

Notice that the condition k ^ 0 assures the existence of the inverse of the 
Hessian matrix. 

The converse of Noether’s theorem, as stated above, shows that all invari¬ 
ants for regular Lagrangian systems are derivable from a Noether symmetry 
which, generally, has non-point character. Hence, for regular Lagrangian 
Ermakov systems, the Ermakov invariant may be viewed as the result of a 
Noether symmetry. In the present case, using the Ermakov invariant (3), 
the Lagrangian (9) and equations (26-27), we find the following Noether 
symmetry associated to the Ermakov invariant /, 

t — t + er, R = R + etR, 9 = 6 + e (— R 2 6/k + t0) , (29) 

where the transformation is expressed in terms of the coordinates ( R,9,t ), 
and r is a completely arbitrary function of time and the dynamical variables. 
The dynamical symmetry (29) gives directly the Ermakov invariant through 
Noether’s theorem. Notice that it is not possible to use the freedom of choice 
of r to reduce (29) to a point symmetry and we necessarily have a dynamical 
symmetry. Consequently the original studies [5]—[8] that considered only 
point symmetries in a one-dimensional Lagrangian formulation, were not 
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apt to explain the Ermakov invariant in terms of Noether symmetries. The 
existence of a (two-dimensional) Lagrangian formalism pointed out in this 
work constitutes an essential ingredient. It is important to note, however, 
that not all Ermakov systems are Lagrangian and, hence the validity of the 
result is still restricted. 

Noether’s theorem may be useful also in the search for second invariants 
for Lagrangian Ermakov systems. With this goal in mind, we restrict the 
treatment to point transformations. In fact, a detailed calculation shows 
that allowing for dynamical symmetries depending at most on velocities does 
not gives additional new results, except for the derivation of the symmetry 
(29). But this symmetry is more easily derived by direct use of the converse 
of Noether’s theorem. However, it is still possible that dynamical Noether 
symmetries depending on higher derivatives may yield new invariants for 
Lagrangian Ermakov systems. 

Let us look for the existence of a second invariant for Lagrangian Ermakov 
systems. For this purpose, Cartesian coordinates (x, y ) are more appropriate 
then ( R,9 ) in the calculation of Noether symmetries. The Lagrangian (9), 
with the corresponding potential (11), becomes 

L = 7j(Ai 2 + IBxi) + Cji 2 ) - (j f(\)d\ +J g(\)d\] - V(R,t ), 

(30) 

where R = R(x,y ) is given by (6). 

Not all Lagrangians of type (30) possess a Noether point symmetry. We, 
therefore, will search for the least restrictive conditions on functions /, g, 
V and on the parameters A, B and C, compatible with the existence of 
Noether point symmetries. In such cases the Noether theorem yields a second 
invariant, independent of the Ermakov invariant. 

The calculation of Noether point symmetries for a given Lagrangian is 
a well known procedure [3] and the main steps of this process can be sum¬ 
marized as follows. Inserting L in the Noether criterion (24) results in a 
polynomial in the velocity components, which must be identically zero. The 
coefficients of the cubic terms imply that 

r = p 2 (f), (31) 

for any arbitrary function p(t) of time only. This information and the coef- 
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ficients of the quadratic terms lead to 


rji = ppx— W(t)(Cy + Bx)+ai(t) , (32) 

V 2 = PPV + W{t) (Ax + By) + a 2 (t) , (33) 

where ai , a 2 and W are arbitrary functions of time. 

From the coefficients of the terms linear in x and y, we get 

dA/dx = Adpi/dt + Bdrj 2 /dt , (34) 

dA/dy = Cdp 2 /dt + Bdpi/dt . (35) 


These last two equations will have a solution A provided that d 2 A/dxdy = 
d 2 A/dydx. This implies, after a simple calculation involving rj\ and 772 given 
by (32-33), that 

W = 0. (36) 

The corresponding solution for (34-35) is 

A = A o(t) + — (pp + P~)R~ T d-i(Ax + By) + a 2 (Bx + Cy) , (37) 


where A 0 (t) depends on time only. The remaining term, independent of 
velocity components in the Noether symmetry condition, yields a first order 
linear partial differential equation for the potential, 


P 


odV OV n dA 

= -2 ppV - 


dt 


+ T) ■ 


dq 


dt 


(38) 


The potential V must comply with the Ermakov form (11), and hence, is not 
completely arbitrary. Inserting ry., r) 2 , A and the Ermakov form for V, we 
get an equation for V, 


2 dV dV 

p w +pp m + 2ppV = 


—A 0 — - (pp + 3 pp)R 2 - di(Ax + By) - a 2 (Bx + Cy)+ 

9 ur ( ry/x px/y \ 

— (ai(Ax + By) + a 2 (Bx + Cy))lJ f(X)d\ + J g(X)dXj + 


( 39 ) 
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K 


( a iy — a 2 x 

V R 2 




Requiring / and g to remain arbitrary implies that the / and g dependent 
parts of the last equation vanish. This is obtained by choosing 


W — ai — a 2 = 0. 


(40) 


Consequently, 

1 

p 2 -fc + PP-Q ^ + 2 ppv = — A 0 - -{p'p + 3pp)R 2 . (41) 

The general solution of (41) is 

= + (42) 

p- Zp p- 

where U(R/p) is an arbitrary function of the indicated argument. In fact, we 
can take A 0 = 0 without any loss of generality since the addition of a function 
of time only in the Lagrangian will not affect the equations of motion. This 
completes the Noether point symmetry calculation. 

We remark that (42) could have been obtained from the potential associ¬ 
ated to the autonomous Hamiltonian Ermakov systems in [14] by performing 
a rescaling transformation x = x/C(t), y = y/C(t), t = / dt/C 2 (t ) and 
requiring the resulting equations to stay autonomous. However, Noether’s 
theorem provides a more consistent basis for the reasons why the potential 
(42) is integrable: integrability here is viewed as resulting from a general 
symmetry consideration and not from an ad hoc transformation. Moreover, 
the calculation of Noether point symmetries provides a hint for the determi¬ 
nation of the more general dynamical Noether symmetries. 

Let us summarize our results so far. From Noether’s symmetry criterion 
applied to the class of Lagrangian Ermakov systems of section 2, we obtained 
the symmetry transformations 

x = x + eppx, y = y-Peppy, i=t + ep 2 . (43) 

The corresponding potentials can be expressed in terms of /, g and another 
arbitrary function U ( R / p) as 

n 1 k ( rv/ x r x /y \ 

v = -^- + - u ( R / f ,) + — U f(\)dX + J g(X)d\J . (44) 
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Consequently four arbitrary functions, namely p(t), U(R/p), f{y/x ) and 
g(x/y), remain in the system. In (R,0) coordinates, the Ermakov invariant 
and the invariant arising from point symmetry have the compact represen¬ 
tation, 

1 . r tan 6 rcot 6 

I = -R^\d)d 2 + J f(X)dX + J g(X)dX , (45) 

J = ^(pR-pR) 2 + U{R/p) + nI{p/R) 2 . (46) 

For the sake of comparison, we writte, using configuration space variables, 
the corresponding non-constant, time-dependent Hamiltonian function which 
follows from (12), 

H =\ te+K 4-' t i + 7 u(R/p) - (47) 

An important point to check is the fact that / and J are in involution. 
This follows from the conventional definition of Poisson brackets and the 
relations x = ( Cp x — Bp y )/n , y = (Ap y — Bp x )/n. 

A simple inspection shows that level surfaces of the function / are not 
compact, a fact that prevents the application of Liouvillc’s theorem. Despite 
this, the exact solution for the Ermakov systems possessing the invariants 
(45-46) can be obtained as follows. For given pit), define new variables 

R = R/p, T=jdt/p 2 (48) 

and transform the Noether invariant (45) into 

J = I ( dR/dT f + U(R) + k I/R 2 , (49) 

an energy-like form. Using the energy method [11], we can find, by quadra¬ 
tures, R as a function of T, and then R as a function of t. The exact solution 
is then obtained by inserting R as a function of t in the Ermakov invari¬ 
ant (45), which can then be interpreted as a first-order ordinary differential 
equation for the solution 6. Being separable, this equation is again soluble 
by quadratures. 

We observe that the Noether symmetry (43) is of the same type as the Lie 
point symmetry eventually admitted by Ermakov systems [23], as expected. 
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Indeed, Noether point symmetries form a subgroup of the more general Lie 
point symmetry group of the same dynamical system [2], At least in this 
case, however, the Noether symmetries result more interesting because they, 
directly, produce constants of motion. 

In spite of the reduction of the Lagrangian Ermakov systems with Noether 
point symmetry to quadrature, these systems are not linearisable, in general. 
Indeed, using (20) and (42) we can show that the only two cases compatible 
with linearisation are 

U(R) = a/R 2 — b/R, (50) 

or 

U{R) = a/R 2 + cR 2 / 2 , (51) 

for a, b and c arbitrary numerical constants. Moreover, when U has the form 

(50) , the function a in (13) must be choosen as a = p. When U has the form 

(51) , a must be choosen as any particular solution for the equation 

« + - p) ^ = ° • (52) 

In all linearisable cases, the function p{t) is left free. 

We finally compare our result with other results in the literature. Gram- 
maticos and Dorizzi have derived the potential (44) with A — C = 1, 
B = 0 and arbitrary p{t) (see equation (60) of [24]) in their search for two- 
dimensional systems with invariants quadratic in velocity. However, they did 
not recognize the Ermakov character of the system. Dhara and Lawande, 
in cases A) and C) in reference [25], found potentials that correspond to 
Ermakov systems or to Noether point symmetry. However, the possible in¬ 
tersection of the two classes of potentials, yielding Ermakov systems with 
Noether point symmetry, was not explored at all. Leach, Lewis and Sarlet 
[26] found a class of potentials with an invariant quadratic in velocity. Within 
this class, there exists an Ermakov subclass, a fact that was not pointed out 
precisely. The Ermakov structure in these systems was hidden due to the non 
usual character of the frequency function, depending on dynamical variables. 
Finally, the second invariant found by Goedert [12] is a Noether invariant 
which can be obtained in the present formalism, by choosing A = C = 0, 
B = 1 and p — 1. 
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4 Conclusion 


In this paper we identified a subclass of Ermakov systems with two inter¬ 
esting properties, namely, of being Lagrangian and of possessing Noether 
symmetries of both point and dynamical character. We also show, by the 
converse of Noether’s theorem, how the Ermakov invariant follows from a 
Noether symmetry of dynamical character. This is an alternative criterion 
to characterize and classify the Ermakov systems which should be consid¬ 
ered parallel with the existence of Lie point symmetries [17, 23] or of special 
structure of the solutions on the complex plane [27]. Moreover, the class of 
Lagrangian Ermakov systems with Noether point symmetries presented in 
this paper may possibly be expanded further by relaxing our starting point 
in terms of the ansatz (9) for the Lagrangian. 
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